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Abstract. We study s-process element abundance ratios in stars by carrying 
out stellar population synthesis, using a rapid synthetic stellar evolution code 
which includes an up-to-date treatment of AGB nucleosynthesis and evolution. 
In contrast to other studies, we find that a large spread in the ^^C efficiency 
parameter (^'^Cos) is not needed to explain the observed spread in the ratios of 
heavy s-process to light s-process elements ([hs/ls]), but this comes naturally from 
the range of different initial stellar masses and their time evolution. As a result, 
the ^^C efficiency needed for fitting most stars in the galactic disk is constrained 
to 1 Ceff < 2. In the same fashion we also study the [Pb/Ce] ratios of lead 
stars and find out that for low metallicities ^'^Ceff ~ 0.5. 



1. Introduction 

5'-cnhanced stars show high abundances 
of elements heavier than Fe compared to 
the Sun, which are produced via slow 
neutron capture processes (s-process el- 
ements). These s-process elements are 
synthesized during the thermally pulsing 
asymptotic giant branch (TP-AGB) phase 
in low and intermediate mass stars, how- 
ever, there are also less evolved stars which 
show over- abundances of these elements, 
thus s-enhanced stars are classified as in- 
trinsic and extrinsic. Intrinsic s-enhanced 
stars (TP-AGB or post-AGB stars) pro- 
duce their own s-process elements includ- 
ing the radioactive element Tc, which is 
observed in their envelopes would have al- 
ready decayed if it was not produced in 
situ. These stars are factories of s-process 
elements in the universe so understanding 



them is crucial to the understanding of 
the chemical evolution of heavy elements 
in galaxies. On the other hand, extrinsic s- 
enhanced stars are not yet evolved enough 
to be in the TP-AGB phase, but have ac- 
creted mass enriched with s-process ele- 
ments from an intrinsic s-enhanced com- 
panion in a binary system. Thus they show 
an overabundance of stable s-process ele- 
ments, but no signatures of Tc. The study 
of these stars is useful for probing the nu- 
cleosynthesis that occurred in their WD 
companions when they were on the TP- 
AGB phase, and it also provides informa- 
tion about stellar interaction in wide bi- 
nary systems, e.g., on the different modes 
of mass transfer and tidal interaction. 

Detailed stellar evolution and nucle- 
osynthesis models must be able to repro- 
duce the abundances of both intrinsic and 
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extrinsic s-enhanccd stars. Testing the ab- 
solute elemental abundances directly turns 
out to be difficult because they depend 
on the dilution of material from the inter- 
shell into the envelope, which in turn de- 
pends on several uncertain features such 
as the amount of dredge up, mass loss, 
mass accretion, etc. However, the s-process 
element abundance ratios are practically 
unaffected by these processes. They pro- 
vide useful constraints on another impor- 
tant uncertain parameter: the total neutron 
flux. This is determined in the models of 
Gallino et al. (1998) by setting the amount 
of ^^C nuclei (^^C efficiency, i.e., ^^Ces), 
which act as the main neutron source via 
the ^^C(a, n)^^0 reaction during the inter- 
pulse periods. GalHno et al. (1998) also in- 
troduce a ^''Ccfr scale where the standard 
case (ST in their paper, i.e., ^'^Cefi = 1) cor- 
responds to the amount of ^^C needed for a 
1.5 Mq AGB star with half solar metallic- 
ity to produce the main s-process element 
abundance of the Sun. It was subsequently 
found that a wide spread in ^^Ceflf (around 
two orders of magnitude) is needed for the 
detailed models to reproduce the spread 
in abundance ratios from the observational 
data (e.g., Busso et al., 2001, Gallino et 
al., 2005, and references therein). However, 
for computing time reasons, these authors 
used only a limited number of initial masses 
for their models and only compared to the 
observations the abundances for a "typi- 
cal" thermal pulse. By carrying out stellar 
population synthesis with a rapid synthetic 
code, we can make models for a much finer 
grid of initial masses and also take into ac- 
count the time evolution of the abundances. 
We can thus produce a much larger set of 
synthetic data, and weight the results with 
the initial mass function (IMF), in order 
to compare to observations in a statistical 
sense. Thus we are able to find better con- 
straints on the ^^C efficiency parameter. 

2. Population synthesis 

We carry out the population synthesis with 
a rapid synthetic evolution code based on 



that from Izzard et al. (2004), but mod- 
ified to be more self-consistent with the 
evolutionary phases prior to the TP-AGB 
(Bonacic et al. in preparation). The s- 
process nucleosynthesis is calculated by in- 
terpolating the results of Gallino et al. 
(1998), using a grid of inter-shell abun- 
dances for different values of stellar mass, 
metallicity, '^'^Ceff and thermal pulse num- 
ber. We investigate a fine grid of initial stel- 
lar masses and we follow the time evolution 
of the surface abundances in each star. We 
explore the ratio of the heavy s-process el- 
ement (hs) average abundance to the light 
s-process element (Is) average abundance 
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in 100 TP-AGB stars for each metallicity in 
a range of initial masses from 1.0 to 8.0 . 
We have taken into account 25 different 
metallicities in the range —2.3 <[Fe/H]< 
0.1 and for each one we have weighted the 
output of the different masses by the IMF 
of Kroupa, Tout & Gilmore (1993). 

The population synthesis is carried out 
for different values of ^^Ceff and the out- 
put is compared to observational data, as 
described below. 

3. Comparison of results with 
observations 

3.1. Intrinsic S-enhanced stars 

In the context of single stellar evolution it 
is straightforward to compare the output of 
the population synthesis with the observed 
intrinsic s-enhanced stars. The first set of 
data that we compare to are the intrinsic 
MS/S stars of Busso et al. (2001), for which 
we select from the synthesized population 
the stars that have a surface effective tem- 
perature less than 3500 K and that are not 
carbon stars, i.e., G/0< 0.95. As can be 
seen in Fig. la, the observations can be fit- 
ted quite well with just one ^^C efficiency 
value, ^^Geff = 2. 
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Fig. 1. Intrinsic s-enhanced stellar population synthesis results compared against the 
observations. The gray scale is a logarithmic measure of the number distribution of 
stars over [hs/ls] or (in panel d) [Zr/Fe]. The crosses are the observational data (see the 
references in the text), which have an average error given by the size of the upper right 
cross in each plot. The number density is weighted by the IMF and by the time each star 
spends in an abundance bin, and then normalized for each metallicity. 



The next set are intrinsic SC and C 
stars, taken from Busso et al. (2001) and 
from Abia et al. (2002). Our synthesized 
SC and C stars are those TP-AGB stars 
which have C/0> 0.95. Again, from Fig. 
lb it can be seen that with only ^''Coff ~ 2 
the data is well fitted, without the need for 
a spread in ^^Ceff. 

The last set of intrinsic stars to com- 
pare to are the post- AGE stars, taken from 
Busso et al. (2001), Van Winckel (2003), 



Giridhar & Arellano (2005) and GaUino 
et al. (2005). We have considered as post- 
AGB stars all those from our synthetic 
TP-AGB sample which have an envelope 
mass of 0.02 Mq or less. Most of them can 
be well fitted with ^-''Ccfr ~ 1 (Fig. Ic), 
with a few exceptions which need ^"^Cefi ~ 
1.5. It is important also to notice that 
there is an apparent split in the observed 
[Zr/Fe], which suggests that some post- 
AGB stars (those with [Zr/Fe] w 0) did not 
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Fig. 2. Population synthesis yield results compared against observations of extrinsic s- 
enhanced stars. The gray scale is a logarithmic distribution of the yield ratios ([hs/ls] and 
[Pb/Ce]) for a population of stars as a function of metallicity. The points are observed 
data from different references (sec text), which have an average error given by the size of 
the upper right cross in each plot. The results are weighted by the IMF and normalized 
for each metallicity. 



experience dredge-up episodes, while oth- 
ers (those highly s-process enhanced) did 
experience dredge- up (Van Winckel, 2003). 
This dichotomy is also observed in our re- 
sults and is consistent with ^^Ceflp « 1 for 
these stars, as can be seen in Fig. Id. 

3.2. Extrinsic S-enhanced stars 

Extrinsic s-cnhanced stars can also be stud- 
ied with single stellar population synthesis 



albeit only in a an approximate way. As 
these stars are formed by accreting mass 
from an s-enhanced companion, their abun- 
dance enhancements are (to first order) 
proportional to the mass of the various el- 
ements yielded by their companions. We 
thus compare to the yields, rather than 
to the surface abundances, of our synthe- 
sized TP- AGB stars, again weighted by the 
IMF. Comparing the population synthesis 
yields with the extrinsic star data of Busso 
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Fig. 3. Range of ^'^Cefi for which the observations of different kinds of s-enhanced stars 
are reproduced by population synthesis and the metallicity range within which they are 
observed. The thick solid lines show the range of ^^Ceflp needed by Gallino et al. (2005) 

to reproduce the observations. 



et al. (2001) and Abia et al. (2002) we 
see in Figs. 2a and 2b that the [hs/ls] ra- 
tios of these stars are well reproduced with 

1 ^ ^■'Ccff ^ 3, which is consistent with the 
values found for the intrinsic s-enhanced 
stars. 

At [Fe/H]< —1 there are no intrinsic 
stars nor is the [hs/ls] ratio very sensi- 
tive to ^'^Ceff changes. It is at this point 
that the role of lead stars becomes impor- 
tant. They arc low metallicity s-enhanced 
extrinsic stars on which Pb has been de- 
tected. Tlic [Pb/hs] ratios are sensitive to 
the value of ^^Cefi-. In particular we com- 
pute the [Pb/Ce] ratios and compare our 
models to the data of Van Eck et al. (2003) 
and references therein. We find that for 
most of the lead star [Pb/Ce] ratios can be 
reproduced with ^^Ces » 0.5 (see Fig. 2c), 
except for two clear outliers. 



4. Conclusions 

We have found that in order to reproduce 

the different [hs/ls] and [Pb/Ce] ratios ob- 
served in s-enhanccd stars only a small 
spread in the ^^C efficiency is needed. The 
observed spread in element ratios can be 
naturally explained by different initial stel- 
lar masses and the time evolution of the 
TP-AGB stars, and it is not necessary to 
use the large spread of ^^Ceff values (more 
than two orders of magnitude) suggested 
by other works, e.g., Busso et al., 2001, 
Gallino et al., 2005, and references therein. 
Galactic disk objects ([Fe/H]> —1) are 
well reproduced by considering the range 
1 ^ ^"^Ccff < 3 in our population syn- 
thesis. Most halo objects (low-metallicity 
stars) are well reproduced with "'^^Ceff ~ 
0.5. This might suggest that at low metal- 
licities ^^Ceff is smaller. Fig. 3 shows a sum- 
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mary of the range of Ccfr values needed to 
reproduce the different types of stellar pop- 
ulation, which is limited within a factor of 
5. 

In the results presented we did not 
take into account the existence of an age- 
metallicity relation, which limits the range 
of stellar masses that actually contributes 
to the observed intrinsic s-enhanced stars 
at each metallicity.In future work we will 
take this into account, and also perform a 
proper binary evolution population synthe- 
sis to get a more reliable outcome for the 
extrinsic s-enhanced stars. 
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